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Abstract 



We present a numerical implementation of the one-loop QED corrections to the hard-bremsstrahlung process e 
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e~e + 7. These corrections can be included in the Monte Carlo event generators employed for simulating Bhabha scattering 
events at low-energy high- luminosity electron-positron colliders. The calculation is performed by employing the reduction 
method developed by Ossola, Papadopoulos and Pittau. Our results are implemented in a modular code for the numerical 
evaluation of the scattering amplitudes for any given phase-space point. In a similar way, we evaluate also the one-loop 
QED corrections to e~e + — > /i~/i + 7, and show an interesting application of the method in the presence of two different 



mass scales in the loops. 
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1. Introduction 

Bhabha scattering (e~e + — > e~e + ) is the process em- 
ployed as a luminosity candle at electron-positron colliders. 
In particular, at high-energy colliders such as LEP, the lu- 
minosity was measured by considering Bhabha scattering 
events at small scattering angles. At high-luminosity me- 
son factories (DA$NE, CESR, etc.) operating at center- 
of-mass energies of 1 — 10 GeV, instead, the luminosity 
is determined by analyzing large-angle Bhabha scattering 
events. In both cases the Bhabha scattering cross section 
is large and dominated by electromagnetic interactions. 
Therefore it is possible to employ the techniques of pertur- 
bative QED and predict the cross section with high accu- 
racy. The smaller theoretical uncertainty which affects the 
resulting cross-section prediction directly translates into a 
more accurate determination of the collider luminosity. 

For phenomenological studies and luminosity determi- 
nations, fixed-order calculations of the Bhabha scatter- 
ing cross section need to be interfaced with sophisticated 
Monte Carlo (MC) generators, which take into account 
realistic experimental cuts and the geometry of the de- 
tectors. The current experimental analyses rely on pro- 
grams such as Babayaga@NLO [H, MCGPJ and 
BHWIDE [21], which include next-to- leading order (NLO) 
corrections and effects related to multiple photon emission. 
The three generators mentioned above agree within 0.1% 
for integrated cross sections and within 1% for distribu- 



tions as shown in Ref. [4[. 

From a fixed-order-calculation perspective, NLO cor- 
rections to e~e + — > e _ e + are well under control. In par- 
ticular, one-loop corrections in the full Standard Model 
were calculated long ago 

NNLO QED corrections to the Bhabha scattering cross 
section play a key role in establishing the accuracy of cur- 
rent MC generators and, eventually, in improving it below 
0.1% as it might be required by electron-positron colliders 
of the next generation, such as the planned International 
Linear Collider. 

It is possible to subdivide NNLO corrections into three 
different sets: i) two-loop corrections to the process e~e+ — > 
e~e + ; ii) one-loop corrections including a single hard pho- 
ton in addition to the outgoing electron-positron pair; Hi) 
tree-level corrections with two hard photons or a hard 
electron-positron pair in addition to the e~e + couple in 
the final state of Bhabha scattering. Obviously, some of 
the diagrams belonging to the sets i) and ii) are infrared 
divergent. However, it is possible in both cases to include 
a class of real soft-photon corrections which can be added 
to the corresponding loop corrections in order to get an 
infrared-finite result for each class. 

The calculation of the two-loop corrections belonging 
to set i) was completed only recently @, 0, H, El- Two- 
photon corrections of class Hi) can be computed with any 
of the existing publicly available tree-level event generators 
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and are safely under control. Radiative corrections of class 
ii) involve pentagon diagrams with one additional hard 
photon in the final state and are not completely known. 
Partial results are available for small-angle Bhabha scat- 
and s— channel annihilation processes at large 



tcring 1C 
angles 

The aim of this letter is to present the calculation 
of the one- loop corrections belonging to set ii). Specif- 
ically, we have realized a modular FORTRAN 95 code 
which allows for a numerical evaluation at fixed phase- 
space points of the one-loop QED corrections to the pro- 
cess e~e + — > e~e + 7 (here 7 denotes a photon with an 
energy larger than a given cut-off threshold) retaining a 
finite electron mass. The calculation has been carried 
out by employing the Ossola-Papadopoulos-Pittau (OPP) 



method 12, 13JI, based on a reduction of the tensor inte- 



grals performed at the integrand level [141 ]. 

As a by-product, we have computed the one-loop QED 
corrections to the annihilation process e~e + — > /i~/i + 7. 
This reaction is an important background for the determi- 
nation of the pion form factor and provides an indepen- 
dent calibration for a measurement of the hadronic produc- 
tion cross section. In addition, it represents an interesting 
application of the OPP reduction method in presence of 
two mass scales inside loop diagrams. Note that a partial 
implementation of NLO radiative corrections to e~e + — ► 
/i~/i + 7 is currently available in PHOKHARA |15f and 
KK MC [II] (see Ref. for a comparison). 

The paper is organized as follows. In Section [2] we 
set our notation and conventions and briefly describe the 
computational technique which we employed. In Section [3] 
we provide numerical results for specific benchmark phase- 
space points. Section[4]contains our conclusions, including 
comments on the numerical stability of the results and on 
the computer time required by the calculation. 



2. Calculation of Radiative Corrections 

2.1. Conventions and Leading Order Cross Section 
We consider the hard-brcmsstrahlung processes 
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fks + fki 



lk 5 



(1) 



where k\ and k 2 denote the momenta of the colliding elec- 
tron and positron, k 3 and k± stand for the momenta of 
the outgoing fermion and anti-fermion, and k$ is the mo- 
mentum of the outgoing hard photon. The particle mo- 
menta obey the mass-shell conditions k\ = k\ 



in. 



^"3 



k\ = m 2 and k\ 



0. The kinematics is described by 
five independent invariants, which can be chosen among 
the six quantities 



(fci + fc 2 ) 2 = AE 2 
{h + h) 2 , 



(2) 
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(ki kj) 2 



(* = 1,2; J = 3,4). 



Here E indicates the beam energy in the center-of-mass 
frame. 








Figure 1: QED tree-level diagrams for e _ e + — > e _ e + 7 (full set) and 
e~ e+ — > /i _ /x+7 (first line). Here and in the following, particles on 
the left side of the diagrams are considered incoming, particles on 
the right side of the diagrams are outgoing. 



The unpolarized leading order (LO) cross section can 
be written as 



(s - Am 2 e ) 
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where Altrcc represents the tree-level amplitude. In the 
context of pure QED, the latter originates from the Feyn- 
man diagrams shown in Fig. [TJ The n-particle phase space 
is defined as 




n+2 



(2^)V 4 > k 1 +k 2 -^2k j .(4) 
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The LO cross section for real-photon emission of Eq. JlJ 



was computed long ago in Refs. [18|, and in the ultra- 
relativistic limit it can be written in a compact way (l9j . 
We report its expression because it can be used for checks, 
although our NLO results, presented in the following sec- 
tions, have been computed without any approximations, 
by retaining the complete dependence on the kinemati- 
cal variables and keeping finite electron and muon masses. 
The ultra-relativistic limit of the LO results [T^| are con- 
veniently expressed in terms of the kinematic invariants of 
Eq. © and of the variables k^ = ki ■ k$, with i = 1, . . . , 4, 

i£|A4 tree | 2 = (4™) 3 T(-^- + 

spins 
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where the function T depends on the final-state fcrmions. 
For the process e~e + — > e~e + 7 one finds 



T = 



1 



ss't 13 t 2 4 



ss'(s 2 + s' 2 ) + t 13 t 24 (t 2 13 + i|t) + 



+£14^23 (^14 + *23) 



while for the case e e + — > [i /i + 7 the function T is 



1 



T — ~i (^13 + *14 + *23 + ^2 
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2.2. Virtual Corrections 

The one-loop corrections to the two scattering pro- 
cesses studied in this work involve a relatively limited 
number of Feynman diagrams: the package QGRAF [20j 
generates 38 one-loop diagrams for the process e~e + — * 



~/i + 7 and 76 diagrams for the process e 



e e + 7. 



Representative graphs are shown in Fig. [21 note that due 
to Furry's theorem, diagrams of class 2c cancel in the sum. 




Figure 2: Representative one-loop diagrams for e e+ 



e e T 7. 



As can be seen from Fig. [3J the calculation of the com- 
plete one-loop corrections requires the evaluation of pen- 
tagon diagrams of class 2g, the most challenging and time 
consuming part of this calculation. 

It is well known that any one-loop amplitude can be 
written as the sum of a linear combination of scalar box, 
triangle, bubble and tadpole integrals, whose analytic ex- 
pressions were obtained in the seminal paper by 't Hooft 
and Veltman [2l|, and a so-called rational part. The ra- 
tional term and the coefficients multiplying the one-loop 
integrals are ratios of polynomials in the masses and the 
kinematical invariants. These coefficients have been tra- 
ditionally determined by means of the Passarino- Veltman 
(PV) reduction method [22j. However, for processes with 
more than four external legs, a straightforward application 
of the PV reduction can generate very large expressions in- 
volving inverse Gram determinants and lead to numerical 
instabilities. 

In recent years, a significant amount of work has been 
devoted to the evaluation of one-loop amplitudes for pro- 
cesses with five or more external legs. Besides improve- 
ments on standard techniques, where a tensor reduction is 
explicitly performed, new unitarity-based numerical and 
analytical developments have emerged (for recent reviews 
see Refs. [23[ where an exhaustive list of references can 
be found). Both directions have been successfully pur- 
sued and have led to impressive results. In particular, 
in the last few months, very challenging calculations in- 
volving four particles in the final state, relevant for LHC 
phenomenology, have been completed 24], 25[. 

In the present letter we describe the evaluation of the 
one-loop corrections to the processes e~e + — > e~e + 7 and 
e~e + — > /i~/i + 7, performed by means of the OPP reduc- 
tion method [12j, ll3j. The one-loop amplitudes, generated 
through QGRAF, have been processed with FORM (2(| 
routines in order to produce a FORTRAN 95 output. 



The latter has then been used as an input for two computer 
implementations of the OPP technique which numerically 
evaluate the coefficients of the basis i nteg rals and the ra- 
tional part for each phase-space point [27J, [28j . Both codes 
have been interfaced with publicly available software which 
allows for the evaluation of the basis integrals 23. 3j 



The one-loop corrections to the cross section can be 
written as 



dRst 



- ]T 2 Re (JWi-ioopJW* ee ) , (8) 



where the interference of the one-loop amplitude A^i-i 00 p 
with the complex-conjugate tree-level amplitude M* ICC is 
evaluated in dimensional regularization and subdivided 
into four contributions, according to the OPP-reduction 
algorithm, 

- 2RC (A-ll-loopXtrcc) =004+11! +K 2 +UV ct .(9) 

spins 

Here CC4 denotes the cut-constructible four-dimensional 
part of the result, written as a linear combinations of scalar 
boxes, triangles, bubbles and tadpoles, IZi + IZ2 stands 
for the so-called rational part and UV c t summarizes all 
contributions induced by the ultraviolet counterterms. As 
a technical remark, we stress that throughout our com- 
putation we have used strictly four-dimensional external 
momenta. 

The rational part has been written through the sum of 
the components 1Z± and 1Z 2 i as explained in Ref. (3lj . The 
integrand corresponding to a generic to— point one-loop 
Feynman diagram can be written as 



N(q) 



a 



' ' 



(10) 



where a bar denotes objects living in d = 4 — 2e space- 
time dimensions, q is the loop momentum, pi are linear 
combinations of the external four-dimensional momenta 
and rrii stand for the masses of the internal legs. The 
numerator N(q) can be split into a four-dimensional and 
an e-dimensional part, denoted in the following by a tilde, 
N(q) = Ni(q) + ^2(5). The four-dimensional numerator 
Ni (q) is then expanded in terms of four-dimensional de- 
nominators Di [l2|: the mismatch in the dimensionality 
between the latter and the d— dimensional inverse propa- 
gators Di of Eq. (TIT))) generates the IZi term in Eq. ©. 
The TZ-2 component, instead, stems from the overlap be- 
tween N2(q) and the ultraviolet poles of one-loop inte- 
grals, and can be evaluated in QED by employing the ad- 
hoc counterterm-like Feynman rules presented in Ref. [3ll ] . 
Note that the overlap between the e-dimensional numera- 
tor N2 (<?) and the infrared poles of one- loop integrals can 
be safely neglected from the very beginning as proven in 
Appendix A of Ref. [H. 

We carried out the ultraviolet renormalization in the 
on-mass-shell scheme: the renormalized charge is chosen 
to be equal to the value of the electromagnetic coupling, as 
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measured in Thomson scattering, at all orders in perturba- 
tion theory; the squared fermion masses are identified with 
the real parts of the poles of the Dyson-resummed propa- 
gators; field-renormalization constants cancel by definition 
external wave-function corrections. 

Throughout our computation we have retained the full 
dependence on the fermion masses, introducing in partic- 
ular the appropriate ultraviolet and 7?-2-type mass coun- 
terterm diagrams depicted in Fig. [31 We observe that, 





Figure 3: Representative mass-counterterm diagrams for e~ e+ — > 
e~ e + 7. Black dots stand for mass-counterterm insertions, necessary 
for performing renormalization and computing the rational term 7^2 
in the massive case. 



in the limit where terms suppressed by positive powers 
of the electron or muon mass are neglected, the contri- 
butions induced by ultraviolet counterterms and 7?.2-like 
rational terms can be embedded in the computation by 
means of simple factors multiplying the tree-level ampli- 
tude. They represent important elements for performing 
checks, as they can be accounted for simply by introducing 
counterterms at the diagrammatic level. Moreover, they 
can be employed for speeding up the numerical evaluation 
of the amplitudes in all cases where terms proportional to 
the fermion masses can be safely neglected. 

By using the U(l) Ward identity in order to eliminate 
the charge counterterm, we find that the ultraviolet coun- 
terterm UV c t in Eq. ^ originates from 



Mi —loop 

[UV 
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47T 



8Z A 



(11) 



where field counterterms arc defined through A p = Z A /2 A^, 

ip = Z^ /2 tp R , Z t = 1 + a/(4ir)5Zi, with i = A,ip, and the 
index R has been introduced for denoting renormalized 
quantities. Explicitly, we have the well-known expressions 
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where the index / runs over the fermions, Nf is the color 
factor (Nf — 1 for leptons and Nf — 3 for quarks), Qf is 
the quantum number associated with the electric charge, 
m/ is the fermion mass, 7 is the Euler-Maschcroni constant 
and fx is the 't Hooft mass. Although QED is not simulta- 
neously free from ultraviolet and infrared divergencies for 
an arbitrary number of space-time dimensions, we have 



followed the common practice of using the same dimen- 
sional regulator for isolating both ultraviolet and infrared 
poles. 

Similarly, within the same approximation, the IZ2 com- 
ponent in Eq. © arises from the interference of A^i_ioo P [^2] 
with the tree-level amplitude, where the former is defined 



as, 



All-loop [K2] = Mtre 
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Y, N fQ 
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(13) 



In both Eq. (|11[) and Eq. (|13[) . vacuum polarization in- 
sertions can be easily removed by switching off all terms 
proportional to the sum over the fermions. 

We observe that the common structure of those terms, 
even in the general kinematic case we worked out, suggests 
that the contributions to the rational term T^i can be fully 
combined with ultraviolet counterterms, thus achieving a 
significant optimization in this part of the calculation. 

The cut-constructible term and the H.\ term in Eq. ([9]) 
have been obtained by direct application of the OPP method. 
Here we describe the two implementations of this part of 
the calculation. 

The first version employs the routines of the publicly 
available package CutTools 27] for the numerical de- 
termination of the coefficients which multiply the basis 
integrals. These routines have been combined with QCD- 
LOOP (2^, used to evaluate the needed scalar integrals 0. 
The basis of loop integrals employed by CutTools in- 
volves also rank-one and rank-two two-point functions in 
order to improve the numerical stability respect to the 
case where the basis contains only scalar integrals, as de- 
scribed in Section 3 of Ref. The additional tensor 
integrals have been included coding the relations provided 
in Ref. 

The second version of the calculation is used for cross- 
checking our results and it makes use of an independent 
FORTRAN 95 code [28] for the reduction of the ten- 
sor integrals, which includes an optional optimization of 
the OPP technique based on the Discrete Fourier Trans- 
form [34[ ■ The basis integrals are evaluated using the code 
OneLOop d(|, written by A. van Hameren. 

It is interesting to observe that the rational part IZi is 
computed by using different strategies in the two versions 
of the calculation: CutTools employs the mass-shift pro- 
cedure introduced in Ref. [Hj], while the second approach 
uses the counterterm-based method described in Ref. [l3| • 

The results obtained with the two independent imple- 
mentations for the cut-constructible term and the compo- 
nent IZi are in very good agreement. 

A final comment concerns the vacuum polarization in- 
sertions (such as diagram 2a in Fig. [2]). Although these 



1 We have used the version 1.0 of CutTools. During the 
completion of this work, version 1.1 has become available at 
http://www.ugr.es/~pittau/CutTools/, with a built-in interface to 
QCDLOOP. 
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types of diagrams are implemented in both codes, they 
are naturally incorporated in any QED computation by 
simply running the fine-structure constant a to the appro- 
priate scale. Furthermore, the introduction of a running 
fine-structure constant allows for a straightforward inclu- 
sion of hadronic contributions using dispersion relations 
and the optical theorem. Since this class of corrections is 
well-known, it will be neglected in the rest of the paper. 

2.3. Real Soft-Photon Corrections 

In order to check our calculation, it is useful to con- 
trol the agreement between the values of the coefficients 
of the residual infrared poles for the virtual one-loop cor- 
rections and those of the analytic expressions of the corre- 
sponding poles arising from real soft-photon emission di- 
agrams. In fact, the sum of the one-loop corrections to 
e~e + — > e~e + 7 (e _ e + — * i-i~fi + "f) and the real-emission 
diagrams for e~e + — ► e~e + j% ott (e~e + — > /i~/i + 7 7 SO ft) is 
infrared finite after integrating over the soft-photon phase 
space up to a given upper cut-off on the energy of the 
undetected photon 7 soft . 

The contribution of real soft-emission diagrams factor- 
izes in the product of the LO cross section of Eq. and 
an infrared-divergent coefficient which can be extracted 
from Ref. fH, 



We define 



acr NLO 



a 



da 



LO 



(14) 



where the functions Jij can be written in the compact form 



- r ^_ ln 1 1 + A. 
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with 



N = 



2 2 

(ki ■ kj) 2 



(15) 



(16) 



and A being a finite remainder which is not considered 
in the present work. Here we have introduced £i = 1 for 
i = 1, 4 and e, = —1 for i = 2,3. In the case i = j, 
obtained by taking the limit /3y — > in Eq. (fT5"|) , one finds 
J ll = 1/2. 

3. Numerical Results 

In this section we show the numerical results for the 
squared LO amplitude, summed and averaged over the 
spins, and the associated one- loop virtual corrections (ex- 
cluding trivial vacuum polarization insertions) for both 
processes at fixed phase-space points. Although the main 
result of our work is represented by a computer code which 
can eventually be interfaced with MC generators, we find 
it useful to show explicitly some numerical results in order 
to allow for detailed comparisons. 



'-LO 



2 

tree | ■> 



(17) 



spins 



(e-%p I* LO = - ^ 2 Re (M 1 - loop Mt rcc ) , (18) 

spins 

where 7 = 0.5772156 • • • and the arbitrary 't Hooft mass 
unit has been set for dcfinitcncss to the value /1 = 1 GeV. 
Following Eq. Q, we further define 

^nlo = ^nlo(cc 4 + n 2 ) +i^ LO (uv ct ), (19) 

isolating the contribution of the ultraviolet counterterms 
from the sum of the cut-constructible and rational parts. 

In the following we will show also the value of the 
residue the infrared pole for real soft-photon emission, in- 
troducing 



Int o — 2lO — ^ ^ Jij , 



NLO 



1 x 

e 



(20) 



with defined in Eq. (Tl"5"|) . 

All results have been obtained using the input data 
suggested by the Particle Data Book [36l | 

3.1. Results for e~e + — ► /i~/i + 7 

We start considering the process e~e + — ► /i~/i + 7. We 
set yfs =10 GeV and show results for the following phase- 
space point: 

fel = (5,0,4.999999973888011,0), 

fc 2 = (5,0,-4.999999973888011,0), 

fc 3 = (0.5738577925797953, 0.4550096995790255, 

-0.3315217520850587,3.476459076282246 ■ 10~ 2 ), 
fc 4 = (4.947897617027284, -1.988746555731907, 

4.244475358996690, -1.581089245425936), 
fc 5 = (4.478244590392921, 1.533736856152882, 

-3.912953606911631, 1.546324654663113), (21) 

where fc, = (k° , kj , kf , kf) and all quantities are expressed 
in GeV. After splitting the one-loop corrections according 
to Eq. (TT9"]) . we obtain: 

Ilo = 5.013964825924999 ■ 10~ 3 , 
I„ LO (CC4 + TCi+TC 2 ) = - 0.3666265876159401 • 10~ 3 + 



+ 1.944055391172138 ■ 10 , 
2"nlo( wv =*) = - - 0.03481919776738505 ■ 10~ 3 + 

-0.3892249340437660 ■ 10~ 3 , 

v 1 t 

I^ LO = - 0.3318073898485551 • 10 3 + 

+ 1.554830457128372 ■ 10~ 3 , 

l5 Tn = - - 0.3318073896945480 ■ 10~ 3 . 



(22) 



All results are expressed in GeV~ 2 and they have been 
obtained by working in double precision. 

In order to check the stability of our results, we have 
performed the so-called N — N test, monitoring the agree- 
ment between the numerical values of the numerator func- 
tion N(q) of Eq. (TUJ)) before and after the decomposition 
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in terms of inverse propagators 27j. When the numerical 
agreement in the comparison does not reach a given limit 
set by the user, the code automatically triggers the use 
of the more time-consuming multi-precision routines 37 1. 



We observe more than 9 digits of agreement between the 
results obtained in double precision, requiring a 10 -5 rel- 
ative precision for the N = N test, and those we got after 
forcing multi-precision in the reduction program for reach- 
ing a 10 -15 relative precision. 

A second test on the precision of our calculation con- 
cerns the cancellation of the poles. After renormalization, 
the residual pole for T^lo ^ s °f P ure infrared origin and it 
matches the infrared pole in 2^(3 w ^ n an agreement of 9 
digits. 

Finally, as stressed in Section FZ72\ we have performed 
the calculation of all contributions by means of two inde- 
pendent codes. 

As a further test of the stability of our results, we study 
the process e~e + — ► /i~/i + 7 for configurations in which the 
muon, or the antimuon, is (almost) parallel to the emitted 
photon. We set y/s = 1 GeV and fix the momenta of the 
colliding leptons and the outgoing photon to be: 

fci = (0.5,0,0,4999997388800458), 
k 2 = (0.5,0,0,-0.4999997388800458), 
fc 5 = (0.4,0.1647604992975971, 

-0.3568677260727233, 0.07415796625373612), 



(23) 



with all momenta given in GeV. 




Figure 4: Virtual corrections for e~e + — » fi~ fi + -f for y/s = 1 GeV as 
a function of the energy of the outgoing muon. The momenta of the 
colliding leptons and the outgoing photon are defined in Eq. (123ft ■ 



Then, we calculate the finite part of the virtual correc- 
tions I^lo °f Eq. (p~8j) , as a function of the energy of the 
outgoing muon £L : the results are illustrated in Fig. [4] 
On the left side of the plot, for low values of £L , we are in 
a configuration in which the muon is almost parallel to the 
photon, while the antimuon goes in the opposite direction; 
we then progressively rotate the muon, keeping the four- 
momentum of the photon fixed, until the muon becomes 
anti-parallel to the photon, while the antimuon becomes 
almost parallel to it. While we observe the appearance of 
collinear divergencies, the virtual corrections I^lo f°ll° w 



a smooth curve, thus suggesting a good stability of our 
result even in this particular kinematic configuration. 



e e + 7 



3.2. Results for e 

Here we consider the process e~e 
yfs = 1 GeV and present numerical results for the follow- 
ing phase-space point: 



e e+7. We set 



fei = 

k 2 = 

k 3 = 

k& — 

k s = 



(0.5, 0, 0.4999997388800458, 0), 
(0.5,0, -0.4999997388800458,0), 
(0.1780937847558600, -0.1279164180985903, 

5.006809884093004 • 10~ 2 , -0.1133477415216646), 
(0.3563944406457374, -2.860530642319879 ■ 10" 2 , 

-0.1832142729949070, -0.3043534176228102), 
(0.4655117745984024, 0.1565217245217891, 

0.1331461741539769, 0.4177011591444748), 



(24) 



with all quantities expressed in GeV. We obtain 

i^ LO (cc 4 .+n 1 +n 2 ) 



0.7586101468103619, 
1 



lNL0("V ct ) 



0.04745064270035045 + 

e 

+0.5005828268263969, 

1 

0.005286348050945757 ■ 

e 

-0.08718044078580632, 

1 

- 0.04216429464940469 + 

6 

+0.4134023860405905, 

1 

0.04216429464958627, 



(25) 



where all results are expressed in GeV~ . 

All the numbers have been obtained by working in dou- 
ble precision and requiring a 10 -5 relative precision for the 
N = N test. Also for this process, we have tested the pre- 
cision of our calculation by forcing multi-precision in the 
reduction program, by checking the complete cancellation 
of ultraviolet and infrared poles and by comparing the two 
independent implementations. 

Following the same idea of Fig. |4l we plot in Fig. [5] 
the finite part of the virtual corrections, I^lo of Eq. (Tl8|) . 
as a function of the energy of the outgoing electron 
ranging between two configurations where the electron mo- 
mentum is (almost) parallel or anti-parallel to the photon 
momentum and collinear divergencies show up. 

4. Conclusions 

We have applied the OPP method to evaluate the com- 
plete NLO virtual QED corrections to the hard brems- 
strahlung processes e~e + — * e~e + 7 and e~e + — * /i~/i + 7, 
both relevant for determining the luminosity at low-energy 
electron-positron colliders. 

It is interesting to note that in the case of muon-pair 
production in association with a hard photon the calcula- 
tion has been performed retaining both the electron and 
muon masses. This shows the capability of the method in 
dealing with several scales inside the loop diagrams. 

The numerical stability of the results has been tested 
by means of three different procedures: the check on the 
quality of the reconstructed numerators given by the N = 



G 




Figure 5: Virtual corrections for e~ e+ — > e _ e+7 for y/s = 1 GeV as 
a function of the energy of the outgoing electron. The momenta of 
the colliding leptons and the outgoing photon are defined in Eq. 1231 1 



N test; the comparison of infrared divergencies arising 
from virtual and real corrections, which cancel as they 
should; the agreement of the results obtained using two 
independent codes. In particular, for LO calculations, we 
find an agreement of at least 12 digits; for NLO results, 
we have 9 digits of accuracy. 

The main result of our calculation is the implemen- 
tation of hard-bremsstrahlung emission corrections, both 
for Bhabha scattering and muon-pair production, into a 
FORTRAN 95 code which employs the publicly available 
package CutTools for the extraction of the coefficients of 
the scalar integrals and the package QCDLOOP for evalu- 
ating the needed scalar integrals. 

The typical order of magnitude of the CPU time of 
the FORTRAN 95 code we have developed is ©(lCT 1 ) 
seconds for each phase-space point on a standard desk- 
top machine. The computational speed, together with the 
modular approach we have followed in computing the ra- 
diative corrections, allows for an implementation of the 
results in the existing MC programs. 

Once supplemented with the finite part for real-emission 
diagrams, the results we have derived allow for a phe- 
nomenological study of the two hard-bremsstrahlung pro- 
cesses e~e + — > e~e + 7 and e~e + — > /i~/i + 7, including the 
calculation of the cross sections and the relevant distribu- 
tions. In addition, they can be compared with the correc- 
tions already implemented in MC generators in order to 
assess the associated theoretical uncertainty. 
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